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Abstract 

For a regular 2n-gon there are (2n —1)!! ways to match and glue the 
2n sides. The Harer-Zagier bivariate generating function enumerates 
the gluings by n and the genus g of the attendant surface and leads to 
a recurrence equation for the counts of gluings with parameters n and 
g. This formula was originally obtained by using the multidimensional 
Gaussian integrals. Soon after Jackson and later Zagier found alter¬ 
native proofs that used the symmetric group characters. In this note 
we give a different, characters-based, proof. Its core is computing and 
marginally inverting Fourier transform of the underlying probability 
measure on S 2 n- Aside from Murnaghan-Nakayama rule for one-hook 
diagrams, the counting techniques we use are of elementary, combina¬ 
torial nature. 


1 Introduction and main results 


Consider a regular, oriented, 2n-gon. There are (2n — 1)!! ways to match and 
glue 2n-sides observing head-to-tail constraint in each glued pair. Each such 
gluing produces an one-face map on an oriented surface. Let eg{n) denote 
the total number of gluings resulting in a surface of genus g. Thirty years 
ago Rarer and Zagier [9] discovered that 
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here n +1 — 25 f is the number of vertices in the map on the surface. Of course, 
eg{n)/{2n — l)\\ is the probability that the uniformly random gluing generates 
a surface of genus g. So, introducing the random variable Gn, n > 1, genus 
of the random surface, so that Vn := u + 1 — Gn is the number of vertices on 
the surface map, and setting Go = 1, we rewrite (Ea as 

/1 I \ y 

1 + 2J2 ^ j . ( 1 . 2 ) 

n=0 ^ ^ 

Their proof used a powerful technique based on multidimensional integrals 
with respect to a Gaussian measure on The identity fll.ip implied a 
remarkable 3-term recurrence for the counts eg{n). A year later Jackson [TT] 
found a group characters-based derivation of an explicit formula for those 
counts, expressed through Stirling numbers of both kinds, and used it to 
prove the recurrence anew. Subsequently Itzykson and Zuber 110] reduced the 
combinatorial calculations in [9]. In 1995 Zagier found another, shorter, proof 
of fll.ll) based on group characters, see [IT] and also Zagier’s Appendix to 
the book [12] by Lando and Zvonkin. In 2001 Lass [13] gave a combinatorial 
derivation based upon the enumeration of arborescences and Euler circuits. 

The Harer-Zagier formula was used by Linial and Nowik [T3] to obtain 
a sharp asymptotic formula for E[G„] and later by Chmutov and Pittel [2] 
to prove that G„ is asymptotically Gaussian with mean {n — logn)/2 and 
variance (logn)/ 2 . Pippenger and Schlech [T3], Gamburd [7], Fleming and 
Pippenger [B] studied the random surface obtained by gluing together edge¬ 
wise n oriented polygonal disks, all with the same number of sides, 3 in [T5] . 
and k > 3 in 0 . i; kn needs to be even, of course. In addition to the 
uniformity of the “gluing” permutation a, Gamburd also assumed that those 
oriented fc-gons were the cycles of the permutation /3 chosen independently 
of a and uniformly at random among all permutations with fc-long cycles 
only, rather than of a hxed such permutation fd. Fleming and Pippenger 
showed that the cyclic structures of the resulting permutations 7 := a/J are 
equidistributed. Gamburd used a Fourier transform-based bound for the 
total variation distance between two probability measures on a hnite group, 
due to Diaconis and Shashahani to prove that, when 2 lcm{ 2 , k] \ kn, 

7 is asymptotically uniform on the alternating subgroup Akn- Fleming and 
Pippenger used Gamburd’s result to obtain very sharp approximations for 
the moments of Vn, conhrming the conjectures made in [15] for the case k = 3. 
Thus the number of vertices Vn in the surface obtained by gluing the given 
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discs and the number of cycles of 7 are equidistributed. Chmutov and Pittel 
[2] extended Gamburd’s result to the general case of n polygons with arbitrary 
“circumferences”, adding up to an even iV —)■ cxd: 7 is asymptotically uniform 
on Ajs! (on resp.) if N—n and N/2 are of the same (opposite resp.) parity. 

In this note we combine Gamburd’s ideas and the Fourier transform on 
S'tv to give a proof of the identity fll.2p . Aside from Murnaghan-Nakayama 
rule, mostly for a simple case of one-hook diagrams, our argument uses only 
elementary enumerative techniques. 

2 Derivation of Rarer-Zagier formula 

Given an even N = 2n, let a and /3 be two independent random permutations 
of [N] chosen uniformly among the permutations with all cycles of length 2, 
and among all {N — 1)! unicyclic permutations respectively. Our task is to 
determine the generating function of the number of cycles of the random 
permutation 7 := a/3. 

The starting point is the Fourier inversion formula for a general proba¬ 
bility measure P on Sn, Diaconis |5]: 

= ( 2 - 1 ) 

■ \\-N 

here A is a generic partition of N, is the irreducible representation of S'tv 
associated with A, = dim(p^), and P{p^) is the matrix-valued Fourier 
transform of P(-) evaluated at p^, 

p(.p") = Y1 pP)PM- 

S^S jy 

Let us evaluate the RHS of fl2.ll) for P = P^, the probability measure on 
Fat induced by a = 11^=1where aj is uniform on a conjugacy class Cj. 
Now P„ = Pa^.^ -k P^^ k ■ ■ ■ k P „^, the convolution of P ^^,..., Pak ■ So, by 
multiplicativity of the Fourier transform for convolutions, 

k 

pap ') = IIGP)- 

j=i 
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Since each is supported by the single conjugacy class Cj, we have (p^) = 
^ Ifx , see [5] . So 


pap "-)=n 4(p")=n 

i=i i=i 

and fl2.1l) becomes 

A j=l 

For a special case s = id fl2.2p becomes 


( 2 . 2 ) 


k 

p(<T=id)= 

A j=l 

and since the LHS is Af{Ci ,..., Ck), the number of ways to write the identity 
permutation as the product of elements of Ci ,..., Ck, divided by 11^=1 l^il; 
we obtain the Fat- version of Frobenius’s identity 

V(C„..., Cj) = )-'■+" n (2-3) 

A j=l 

Zagier’s proof in im used (12.dh for fc = 3. In our argument we use (12.21) for 
k = 2 only, heavily relying instead on arbitrariness of s G S^. For ai = a, 
(J 2 = (3 and a = ''y = af3, this equation becomes 

M (2-4) 

A 

C 2 {Cn resp.) consists of all {N — 1)!! permutations with cycles of length 2 
only (all {N — 1)! unicyclic permutations resp.). 

Let Ci^ be a generic conjugacy class comprising all permutations with 
given counts Uj of cycles of length j G Dehne J = Jiv) = {] G 
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[l,N] : Uj > 0}. Let ck = {ai,a 2 , ■ ■ ■) be an arbitrary composition of N 
containing Vj components eqnal j, {j G J). By Mnrnaghan-Nakayama rnle, 
Stanley [16] (Section 7.17, Eqnation (7.75)), 

P(C) = (2.5) 

T 


where the snm is over all rim hook diagrams T of shape A and type ex, i. 
e. over all ways to empty the diagram A by snccessive deletions of the rim 
hooks, one hook at a time, of lengths ai, a ^,.... Fnrther ht(T) is the snm 
of heights of the individnal hooks (nnmber of rows minns 1) in the hook 
diagram T. 

This remarkable formnla implies that ~x^{Cfq') = 0 unless the diagram A 
is a single hook A*, with one row of length Ai and one column of height A^, 
Ai + A^ = + 1, in which case 

X\C^) = (-1)"'-'. (2.6) 


Consider x^*(C 2 ). If Ai is even, then A^ is odd, so that for each T the rim 
hook deleted last consists of two hrst cells in the row. So ht(T) = (A^ — l)/2 
for every T, and the number of T’s is the number of ways to intersperse 
first (Ai — 2)/2 deletions, from left to right, of domino tiles in the row with 
(A^ — l)/2 deletions of domino tiles, from bottom to top, in the column, and 
this number is 


/ Aj—2+A^ —1 \ / N—2 \ / N—2 \ 


■ N^2 
2 

Ai-2 

2 


Analogously, if Ai is odd then 

/ N-2 ' 

Combining fl2.6p - fl2.8p . we have 

.y(C2)x"‘(c«) = (-i)‘"‘-'>'"(^) 

/ N-2 \ 

(C 2 )x^‘ (Cw) = (-l)C'+2)/2 j , if Ai is odd. 


if A 1 is even; 


(2.8) 


(2.9) 
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As for /^*, applying the Hook formula we obtain 




N\ 




N-1 

Ai-1 


Introducing 


F(1V.A.):= 

we transform fl2.ip into 


if A 1 is even; 


if Ai is odd, 


= WAi)x'''(s); 

^ ' Ai = l 


( 2 . 11 ) 


( 2 . 12 ) 


here A* is a hook of size N, with the row of length Ai, and the column of 
height A^ = + 1 — Ai- 

Turn to X^* i^)- Let s G Ci^. Introduce i = l{u) := min J = min{j : i/j > 
0 }, the size of the smallest positive component of a. 

Case Ai < Here, for every ck, there is only one rim hook tableau T, 
since the arm of the hook deleted last must be the whole row of A*, and, 
denoting v = ^j', 

ht(T) = \^-u = {\^ + u) (mod 2) ^ x^\s) = (-1)^'+A (2.13) 


Case Ai > Since in fl2.5p the composition a can be chosen arbitrarily, 
let us assume that components of o: are non-increasing. So the composition ol 
consists of the segment formed by all components of the largest size, followed 
by the segment formed by all components of the second largest size, and so on, 
all the way to the terminal segment formed by all components of the smallest 
size £. Therefore in a generic tableau T the last rim hook is a diagram /i* 
comprising one row and one column of sizes fxi and /i', pi -|- /i^ = £ + 1. 
All the other rim hooks in T are either horizontal or vertical, successively 
deleted from the row and from the column of A* respectively. Let hr be the 
number of those horizontal hooks of size r G J; so G [0, i>r] for r > £, and 
hi G [0, z/£ —1]. The admissible h = {hr} must meet the additional constraint 

'^Kr + fii = Xi- (2.14) 

reJ 
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The total number of the tableaux T, with parameters /ii, h, is 



Further, the p^-long leg of the last hook contributes /i^ — 1 to the height 
ht(T), while the total contribution to ht(T) of those vertical rim hooks is 
(A^—the sum of their sizes, minus — hr), their total number. 

So 


ht(T) = (/i^ - 1) + (A^ - /i') - 5^ {l^r - - hr) 

reJ 

= A^ — 1 + {ur — hi^r — hr) (mod 2). 

rej 


Therefore the total contribution to (s) of the rim hook tableaux T with 
the last rim hook fi* is 


^ Y^(^_iyr-Se,r-hr I ^i^r 

h meets 112.1411 rgj 


hr 


'y ' • ^ ^_yur-Si^r-hr 
r^J hr 


— hl^r 

hr 


r£j 


To get x^*('S) we need to sum this expression for all 1 < /xi < min{f', Ai}, 
i. e. for 1 < /ii < £, because Ai > 1. For those /xi, Ai —/xi ranges from Ai — £ 
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to Ai — 1. Thus 


Ai-l 


x"'W = (-i)"‘-‘ K'l nK'-1)’"'*" 

t=Xi—i rGJ 




EnllK'-i)' 

r=0 / r£j 

!!(?’■ - 1 )"-'* 

r&J 




r£j 


,_i)A.+.[jA.)^n(i-e) 


(2.16) 


rdJ 


Observe that, for Ai < £, we have r > Ai — 1 for all r G J; so the bottom 
RHS in fl2.15p is 


(_l)A.+.|^A.-.| (1 _ 

reJ 


which is the value of {^) see (I2.1,3p . Therefore, the bottom line 

expression in (l 2 .L^p for holds for all Ai. 

Though seemingly unwieldy, the formula fl2.12p for P (7 = s) together 
with fl2.15p instantly lead to a promising, intermediate, expression for the 
marginal distribution of X„, the number of cycles of 7 . 

Theorem 2.1. For i/ G [1, iV], iV = 2n, 


P{X^ = u) 


N 


Ai=l 


X 


1-e vi-^e 


(2.16) 


or equivalently 


N 


£■[9-^"] = E(-Tm.Ai) • 

Ai=l 






y-\ 


(2.17) 






















Proof. By fl2.12p . 


P(X„ = u) 


1 ^ 

Ai = l s: u{s)=i' 


Here, by fl2.15p and Cauchy formula \C^\ = N\ , , we obtain 


E x"‘w = (-i)"‘iv!r‘i^ n 

s:iy(s)=iy Ylr^r=t^, r>l 

rur = N 


{i-OYl 


Url 


= (-i)^‘iv!|y>]^[iVin 

^ r>l 


E 

Ur>(} 

„-]exp[-J] 

\ r>l 


■yx^{l-n 


yx^{l-C) 


A 


= (-l)^‘Af! ly’lY"^ [x'^y"] exp f-y log + »log ^ 


= (-i)^‘Af!K^>lY4j[^V] 


1 — X 
1 


The proof is complete. □ 

Let us use (I2.17p to get an explicit formula for P(X„ = i'). To begin, 
from fl2.1ip it follows that, both for Ai = 2m, (0 < m < u) and Ai = 2m +1, 
(0 < m < n — 1), 


(-i)"‘mA,) = (-i) 


n—m+1 




Q{n,m) : = 


\2m-l) 

(2m — 1)!! (2(n — m) — l)!! 

( 2 n- 1 )!! ^ 


(2.18) 


with (—1)!! := 1. Next we evaluate H(n, z/, Ai), the coefficient of 

in the expansion of (A^) • known, Comtet [3] (Section 5.5), 

that 


1 

V\ 



= A 


s(£, z/) 

p. 


(2.19) 


9 





























where s{i, v) is the (first kind Stirling) number of permutations of [£] with v 
cycles. So, setting 


1 — X 1 — 1) 

1 — 1 — rj ^ 1 — X ' 


we have 



Next, for £ < 2n, 


\x 


^2ni 


1 ) 

1 — X 


= \x 


2n-l-\ 


e-1 


X 


= {i-1 )^ 


2 n-£ 


(- 1 ) 


2n-l 


and finally 

- 1 )'=- !)'-■ = ■ 


Collecting the pieces we arrive at 


A(n, z/, Ai) 



2 n-l\ ^ s{l,v) [2n-\-\ 

^ ^ l\ \2n-t)- 


( 2 . 20 ) 


Let Pe{Xn = I') and Po(Xn = i') denote, respectively, the contribution of 
even Ai and odd Ai to the RHS of fl2.16l) . To compute Pe{Xn = v) and 
Po{Xn = i') we will need two simple identities 





'\2a—h 


(- 1 )^ 


h — a 


6 _|_ 


6 — a — 1 




-h 
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directly implied by 


6>0 j=0 


6>0 j=0 


'a\ [ 2 j' 




= [!-(! + = (-l)“x“(2 + x)“, 



= (-l)“(a;“ + a;“+')(2 + x)E 


(see Gould [8] or www.matli.wvu.edu/gould/, Vol.d.PDF, (1.62)). By fl2.20p . 
and fl2.18p . 


Fe(X„ = z/) = 5^(-l)’^-”*g(n, 

t=i 

E(-i) 


m] 


m=l 

X 


2n- 1 


2 (n — m) 

^ «(£, z/) (2{n — m) 


t>u 


l\ 

s(£, v) 


2n — 


E(-i)'^E(-i) 


1>U 


3=0 


n — 1 


2 j 
2n — 


E(-i) 


v) (n- I 


l>u 


£! V £ - 2 


Similarly 


£> 1 / 


i=o 


j J \2n — ij 


E(-i) 


e>u 


e+n-l Sj^: 

a 


.£-2 


n — 1 
i-2 


+ 2 


£-1 


Combining fl2.2ip and fl2.22p we have 

\£+n-lo£-l'®(^’ 


P(X„ = Z/) = ^(-1)^+-12' 


1>U 


e. 


n — 1 
-2 


+ 


rz, — 1 


n — 1 
l-l 


^£+n-l9£-l^(^G ^^ 

’ £! U -1 

£>y ^ 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 
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(This simple formula appears to be new.) Therefore 



1>V 


1 ^ 1 




2 


2 


by fl2.19p . as ^ = 1 — Consequently 



whence, setting Xq = 1 , 



which is the Harer-Zagier formula for the genus as X„ = n + 1 — 2G„. 

Acknowledgment. I owe a debt of gratitude to Huseyin Acan for innu¬ 
merable discussions of the random chord diagrams. It is my pleasure to 
thank Sergei Chmutov and Sergei Duzhin for introducing me to the remark¬ 
able Harer-Zagier formula. I am grateful to Sergei Chmutov for patiently 
explaining to me the topological issues of the gluing models, and for a very 
helpful feedback on an initial draft of this paper. 

References 

[1] S. Chmutov and B. Pittel, The genus of a random chord diagram is 
asymptotically normal, J. Combin. Theory Ser. A, 120(1) (2013) 102- 


110 . 


[2] S. Chmutov and B. Pittel, On a surface formed by randomly gluing 
together polygonal discs, Adv. Appl. Math, (submitted). 

[3] L. Comtet, Advanced Combinatorics, Dordrecht-Holland, (1974). 

[4] P. Diaconis and M. Shashahani, Generating a random permutation with 
random transpositions, Z. Wahr. Verw. Gebiete, 57 (1981) 159-179. 


12 
















[5] P. Diaconis, Group Representations in Probability and Statistics, IMS 
Lecture Notes-Monograph Series, 11 (1988). 

[6] K. Fleming and N. Pippenger, Large deviations and moments for the 
Euler characteristic of a random surface, Random Struct. Algorithms, 
37 (2010) 465-476. 

[7] A. Gamburd, Poisson-Dirichlet distribution for random Belyi surfaees. 
The Ann. Probability., 34 (2006) 1827-1848. 

[8] H. W. Gould, Combinatorial identities; a standardized set of tables list¬ 
ing 500 binomial coefficient summations, 2nd Edition, Morgantown, W. 
Va. (1972). 

[9] J. Rarer and D. Zagier, The Euler characteristic of the moduli space of 
curves. Invent. Math., 85 (1986) 457-485. 

[10] G. Itzykson and J. B. Zuber, Matrix integration and combinatorics of 
modular groups, Gomm. Math. Phys. 134 (1990) 197-207. 

[11] D. M. Jackson, Counting cycles in permutations by group characters, 
with an application to a topological problem, Trans. AMS, 299 (1987) 
785-801. 

[12] S. K. Lando and A. K. Zvonkin, Graphs on Surfaces and Their Applica¬ 
tions, (Encyclopaedia of Mathematical Sciences, volume 141), Springer- 
Verlag (2004). 

[13] B. Lass, Demonstration combinatoire de la formule de Harer-Zagier, G. 
R. Acad. Sci. Paris, 333 (2001) 155-160. 

[14] N. Linial and T. Nowik, The expected genus of a random chord diagram. 
Discrete Gomput. Geom., 45 (2011) 161-180. 

[15] N. Pippenger and K. Schleich, Topological characteristics of random tri¬ 
angulated surfaces. Random Struct. Algorithms, 28 (2006) 247-288. 

[16] R. P. Stanley, Enumerative Combinatorics, Vol. 2, Gambridge Univer¬ 
sity Press (1999). 

[17] D. Zagier, On the distribution of the number of cycles of elements in 
symmetric groups, Nieuw Arch. Wiskd., 13 (1995) 489-495. 


13 



